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Abstract: We investigate a connection between generalized Fibonacci 
numbers and renewal theory for stochastic processes. Using Blackwell's re- 
newal theorem we find an approximation to the generalized Fibonacci num- 
bers. With the help of error estimates in the renewal theorem we figure out 
an explicit representation. 

Keywords: generalized Fibonacci numbers, renewal theory 
Mathematics Subject Classification (2000): 11B39, 60K20 

1 Introduction 

A variety of generalizations for the Fibonacci numbers has been proposed 
over the last decades and were studied using different techniques. In this 
article we study the (^-generalized Fibonacci numbers (d > 2) defined as 
follows: 

= for all n < 0, F[ d) = 1 

and by the recursion 

d 

^ = E^2i for all n> 2. (1) 

Obviously the standard Fibonacci numbers come out for d = 2. The case 
d = 3 is also of special interest, the resulting sequence is called Tribonacci 
sequence and the first few terms for n > are given by 

0, 1, 1, 2, 4, 7, 13, 24, 44, 81, 149,... 

The generalized Fibonacci numbers were introduced by E. P. Miles in [5] and 
were further studied e.g. in [3] and [4], where Binet-type representations were 
given. These kind of representations were obtained using different standard 
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techniques for difference equations such as generating functions and matrix 
techniques. From these techniques one expects to obtain representations of 
the form 

d 

F<f> = for all n > 0, (2) 

i=i 

where z%, z& are the (possibly complex) roots of the polynomial x d — x d ~ l — 
... — x — 1. The coefficients Ai, A^ are not so straightforward to find. 

In this article we take a stochastic point of view to represent the general- 
ized Fibonacci numbers via a stochastic process and use Blackwell's renewal 
theorem we find the asymptotic behavior of the generalized Fibonacci num- 
bers using one special root of the above polynomial in Section [2j Here it is 
remarkable that the right constant comes out immediately. 
Since we wish to get an explicit representation of for all n we use geo- 
metric error estimates for the renewal theorem to see that our approximation 
is sufficiently good in Section [3j In opposite to the representation in (jSJ) our 
representations is just based on one special root of the polynomial. 



2 Approximation using renewal theory 

First we give a combinatorial interpretation of the (i-generalized Fibonacci 
numbers. Fn is the numbers of possibilities to write n — 1 as the sum of 
numbers in {1, d}\ 

Lemma 2.1. With the notation 

oo m 
m=0 k=l 

it holds that 

FW = \jM 1 \ forallneZ. 

Proof. For all n < it holds that Fn = = |0| = (J 7 ^!! and for n = 1 we 
have F^ = 1 = |{()}| = |^o^|- For n > 1 the following recursion holds 

d m 

l-^n^l = I [J{( x i)i==l ^ '■ x m = k}\ = l^n-fcl' 

k=l k=l 

which proves the claim. □ 
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For the use of renewal theory we embed the generalized Fibonacci num- 
bers into a stochastic process as follows: 

Let Xx,X 2 ,... be independent and identically distributed random variables 
with P(X\ = i) = q l for all i = 1, d, where q = qd is the unique number 
in (0, 1) with 

q + q 2 + ...+q d = 1 (3) 

and let Sk = J2i=i Xi,k = 1, 2, ... be the random walk generated by X 1 , X 2 , .... 
Furthermore let r n = inf{k > : Sk > n} denote the first passage time over 
n for n = 1,2, .... 

Remark 2.2. Formally for d — oo the random variables Xi,X 2 ,... have 
a geometric distribution with parameter 1/2. For geometric distributions 
renewal theory turns out to be much easier than in the general case. 

Remark 2.3. By multiplying equation (TJ|) by l/q d we see that 1/q is a root 
of the polynomial x d — x — ... — x — 1. 



Proposition 2.4. For all n > 1 it holds that 

P(S Tn =n)= q n F { n %. 
Proof. For all ...,x m ) G jffl we have 



ppri = X1 , ...,x m = Xm ) = n p (^ = ^ = 



i=l 



hence using Lemma I2TT1 



□ 



The proposition above states that we can represent in terms of the 
probability that Si, S 2 , ... visits the state n. This probability is well studied in 
probability theory, more specifically in renewal theory. The random variables 
Xi,X 2 ,X 3 , ... are interpreted as renewal lifetimes and P(S Tn = n) can be 
seen as the probability for a renewal at time n. The asymptotic behavior of 
P(S Tn = n) for n — > oo is well understood: 

The Blackwell renewal theorem for the lattice case (cf. e.g. [U Theorem 
1.2.2]) states that 

P(S n) : 1 1 (9 ' 1)2 

for n — > oo. This leads to the following 
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Theorem 2.5. 

F {d) ^ Crfg -(n-l) ) 

where ~ means asymptotic equivalence. 

Proof. Using Lemma [2.11 it holds that F^ 1 q n = P(S Tn = n) — > c&. □ 

Remark 2.6. With equation (TJ|) m mine? it is not /iard to believe that 
converges to a constant, but it is not clear how to find the right constant eas- 
ily. As seen above for our approach nearly no any algebra is needed. 

3 Error estimation 

When using the approximation c d q'^ n ~^ for it is first of all not clear 
if the approximation is good or not. As an example let us consider the 
Tribonacci sequence 

0, 1, 1, 2, 4, 7, 13, 24, 44, 81, 149,... 

In this case C3g _( - n_1 ) is given by 

0.33, 0.61, 1.13, 2.09, 3.84, 7.07, 13.01, 23.94, 44.03, 80.99, 148.98,... 

and we see that the approximation is close to the true value; the error is 
always less than 1/2, so that rounding the approximation gives the correct 
value. 

In this section we estimate the error in the general case when using the 
approximation to the generalized Fibonacci numbers obtained in the previous 
section, i.e. we study 

x n = F«> - c d q- [n - l) = g _(n - 1) [i 3 (5' Tn _ 1 = n - 1) - c d ], n E Z. (4) 

Because we want to find a global bound for \x n \ we especially have to show 
that \P(S Tn = n) — a\ has a geometric convergence rate. The question of 
convergence speed in the renewal theorem has been studied intensively and 
we apply these results in our setting. To this end we first show that X\ is 
"new better then used": 

Lemma 3.1. For all i,j > with P{X\ > i) > it holds that 
P(X 1 >i+j\X 1 >i)<P(X 1 >j). 
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Proof. The proof is just straightforward calculus: 

i i+j 



P(X 1 >t + j) = l-P(X 1 <t+j) = l-J2Q k - Yl 9* 



k=l k=i+l 

and 



fc=i+l fe=l k=l 



P{X X > i)P{X 1 > j) = (l-£ g fe )(l-^ g fe ) = 9*) E 

fe=l fc=i fc=l fc=l k=l 

therefore P{X\ > i + j\Xi > i) < P(X\ > j) is equivalent to 

i+j j i 

i.e. equivalent to 

> i - J2 qk = J2 qk = qt J2 q 

k=l k=i+l k=l 

what is obviously true. □ 

Now we can prove the following 
Theorem 3.2. 

/ 1 \ n—\ 

\x n \ < (1 — q) I I for all n > 1. 



k 



Proof. For the proof we use the results from [2] on the rate of convergence in 
the renewal theorem: Using Lemma [3. II we can apply [2, Corollary 2.1] and 
obtain 

{PiSr^ = n - 1) - c d \ < P{Xx > l) n for all n > 1, 



hence 



n | = g -(»- 1 )|P(5' TB _ 1 =7i-l)-c| 

< q- {n - 1) P(X 1 > l) n = g-( n - x )(l - g) n 



1 - g 



for each n > 1. □ 
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As a result we obtain the representation of Fn in terms of Cd and q as 
follows: 

Corollary 3.3. \x n \ < 1/2 for all n G Z, so that is the unique integer 
with minimal distance to c^g"" 1 . 

Proof. Since q > 1/2 it holds that (1 — q)/q < 1, so by (i) we have < 1/2 
for all n > 1. Hence it remains to show \x n \ < 1/2 for all n < 0. Since in 
this case |x n | = caq - ^ 1 ' we only have to treat |xo| = q/E(Xi). We have 

d 

= > l + 2q + q 2 + ... + q d - 1 

> l + g + g 2 + ... + g d - 1 + g d = 2, 
so that \xq\ < 1/2. □ 

Remark 3.4. We just want to mention that the previous result can also 
be obtained elementary without the use of renewal theory. Furthermore the 
representation given above is different to the usual Binet-type representation 
since it just uses one root of the characteristic polynomial. This has the 
benefit that not all (complex) roots of the polynomial have to be found. 

Remark 3.5. With equation (TJ|) in mind we can furthermore see that 1/q is 
the only root of the equation x d — x d ~ x — ... — x — 1 with modulus > 1, because 
else x n = F^ — Cdq n ~ l could not converge to 0. This is one more proof for 
this well known result. 
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